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1 .  Abstract 


In  this  project,  we  have  focused  on  how  to  improve  the  thermoelectric  properties  of  Silion  nanowires,  in 
particular  by  using  doping.  We  have  studied  two  kinds  of  doping  for  Silicon  nanowires: 


(1 )  Silicon  nanowire  is  doped  with  its  isotope.  We  found  that  when  ^®Si  nanowire  is  doped  with  50%  ^®Si, 
its  figure  of  merit  (ZT)  can  be  enhanced  by  31%.  The  results  have  been  published  in  Applied  Physics 
Letter  95,  063102  (2009); 


(2) 


Silicon  nanowire  is  doped  with  Ge.  Our  study  shows  that  when  n-type  Siliocon  nanowire  is  doped  with 
50%  Ge,  its  ZT  can  be  4.3  times  of  that  of  Si  nanowires.  This  means  that  if  the  ZT  value  of  Si  nanowire 


is  about  0.6-1  as  claimed  by  experiments,  then  the  ZT  value  of  n-type 


can  be  as  high  as 


2. 5-4. 3.  This  makes  Ge  doped  Si  nanowire  a  promising  candidate  for  industrial  applications.  The 
results  have  been  submitted  to  Applied  Physics  Letter  for  publication. 


2.  Objectives: 

The  objective  of  this  proposal  is  to  study  the  possible  ways  to  improve  ZT  of  semiconductor  nanowires 
by  reducing  thermal  conductivity  from  phonons.  We  will  concentrate  on  the  computational  study  of  how 
isotope  (impurity),  quasi-periodic  structure  (interface)  and  the  substrate  reduce  the  thermal  conductivity 
due  to  phonons  thus  increase  the  ZT  of  different  semiconductor  materials. 

3.  Introduction: 


Thermoelectric  material  will  play  an  important  role  in  recycling  the  wasted  heat.  The  efficiency  of 
thermoelectric  material  is  determined  by  the  figure  of  merit:  vj-,  here  S  is  the  Seebeck  coefficient, 

K 

a  is  the  electrical  conductivity,  7  is  the  absolute  temperature,  and  k  is  the  thermal  conductivity,  k  =  Ke+  Kp, 
where  Kg  and  Kp  are  the  electron  and  phonon  (lattice  vibration)  contribution  to  the  thermal  conductivity, 
respectively.  In  order  to  make  thermoelectric  materials  competitive,  the  ZT  of  the  material  must  be  larger 
than  three.  However,  the  road  to  achieve  this  goal  has  not  been  very  smooth  so  far.  First,  simple  increase 
S  for  general  materials  will  lead  to  a  simultaneous  decrease  in  a.  Also,  an  increase  in  a  leads  to  a 
comparable  increase  in  the  electronic  contribution  to  k. 

Nanostructured  materials  have  been  come  to  the  spotlight  in  the  last  decade.  It  is  expected  to  have 
larger  ZT  compared  with  the  bulk  counterpart.  The  possible  underlying  mechanism  is  that  the  density  of 
states  near  the  Fermi  energy  level  is  increased  and  the  reduction  in  the  thermal  conductivity  can  be 
reduced  significantly  [1-4].  For  example,  Silicon  nanowire  (Si  NW)’s  thermal  conductivity  is  foundlOO  times 
smaller  than  that  of  bulk  silicon  [5-10]  because  of  the  rough  surface  scattering. 

Moreover,  our  molecular  dynamics  (MD)  modelling  demonstrate  that  isotopic  doping  provides  an 
effective  approach  to  decrease  thermal  conductivity  of  carbon  nanotube  by  50%  [11],  this  theoretical 
prediction  has  been  confirmed  experimentally  by  Berkeley’s  group  [12].  Therefore,  we  can  also  apply  the 
same  strategy  to  reduce  further  the  NW’s  thermal  conductivity.  Indeed,  in  our  recent  simulation  [13],  we 
found  that  the  thermal  conductivity  can  be  tuned  by  the  Ge  contents  thus  open  the  door  for  increasing  the 
ZT  further.  Indeed,  it  has  been  observed  experimentally  that  when  the  Ge  content  in 

nanocomposites  increases  from  5%  to  20%,  the  thermal  conductivity  decreases  obviously.  [14]  However, 
at  the  same  time,  the  power  factor  also  decreases,  and  induces  uncertainty  in  the  change  of  figure  of  merit 
ZT.  Although  a  number  of  studies  have  been  reported  recently  about  the  thermal  and  thermoelectric 
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properties  of  nanocomposites  and  nanowire[13-20],  it  is  still  not  clear  on  the  quantitative  impact  of 

the  composition  on  the  thermoelectric  performance. 

In  this  project,  we  would  like  to  investigate  the  doping  effect  Germanium  to  Silicon  nanowire. 

Since  the  big  difference  of  the  Germanium  atoms  from  Silicon  atoms,  we  should  expect  a  large  reduction  of 
thermal  conductivity  of  Siliocn  nanowire  more  than  50%,  thus  can  double  the  ZT  of  silicon  nanowire.  More 
precisely,  we  would  like  to  calculate  the  doing  effect  on  all  relevant  physical  quantities  like  Seebeck 
coefficient,  thermal  power,  the  electrical  conductivity  and  thermal  conductivity  due  to  both  electrons  and 
phonons,  and  eventually  ZT .  With  computer  simulation,  we  expect  better  understand  of  the  underlying 
physics  thus  we  shall  be  able  to  work  out  a  recipe  for  the  design  of  the  silicon  nanowire  for  the 
thermoelectric  application. 


4.  Theory  and  Computational  methods 


The  electrical  conductivity  o,  the  thermal  conductivity  due  to  electrons  Xg,  and  the  Seebeck  coefficient  S, 
are  obtained  from  the  electronic  structure  with  the  solution  of  1 -Dimensional  Boltzmann  transport  equation 
as: 


e  T 

5  =  J_(A<“^)-'A<'^ 


^(«) 


Yaa 


/?exp(;g(A,  -/u)) 

(l  +  exp(/?(A^-//)))^ 


p{E,)E,{E,-My 


Here  e  is  the  charge  of  carriers,  T  is  the  temperature,  is  the  electron  energy,  t  is  the  relaxation  time,  m* 
is  the  effective  mass  of  the  charge  carrier,  p  is  the  electron  chemical  potential  and  D(Ek)  is  the  density  of 
states.  The  relaxation  time  t,  is  a  complex  function  of  the  electronic  structure,  temperature  and  carrier 
concentrations.  The  purpose  of  our  work  is  to  determine  the  qualitative  effect  of  SiNW’s  transverse 
dimension  on  its  thermoelectric  property.  In  this  work,  we  obtain  the  value  of  relaxation  time  by  fitting  the 
calculated  mobility  to  measured  electrical  conductivity  data  of  SiNW. 

As  the  lack  of  experimental  electrical  conductivity  data  on  diameter  dependence,  here  we  use  the 
experimental  data  of  SiNW  with  fixed  diameter  of  48  nm  (n=1 .7x10^®  cm‘^,  a=588  (Qcm)  \  from  Ref.  [5], 

and  obtain  r  =  4.3x10“''’ s  for  n=1.7xl0^®  cm‘^.  Then,  we  use  the  dependence  relation  between  the 
mobility  and  carrier  concentration  in  bulk  silicon  in  Ref.  [21]  to  calculate  the  carrier  concentration 
dependent  relaxation  time.  For  example,  r  =  1.7x10“'^  s  for  n=1. 0x10^^  cm‘^,  and  r  =  9.7x10“'®  for 
n=5.0x10^^  cm'^.  As  ZT  of  n-doped  SiNWs  is  considerably  larger  than  that  of  their  p-doped  counterparts, 
we  only  study  n-doped  wires.  The  carrier  concentration  is  defined  as: «  =  | D(E  - //) x  f{E - /j)xdE ,  where 

f(E)  is  the  Fermi  distribution  function. 


For  S^_^Ge^N\Ns  we  focus  on  that  NWs  oriented  along  the  [110]  direction  with  rectangular  cross 

section  area  of  2.3  nm^.  The  atomic  structure  is  initially  constructed  from  diamond  structured  bulk  silicon. 
Then  Si  atoms  are  randomly  substituted  by  Ge  atoms  and  the  geometry  is  relaxed  to  its  closest  minimum 
total  energy.  The  surface  dangling  bonds  are  terminated  with  hydrogen  atoms.  A  supercell  approach  is 
adopted  where  each  wire  is  periodically  repeated  along  the  growth  direction.  The  size  of  the  supercells  in 
the  transverse  plane  is  large  enough  (>15  A  from  surface  to  surface).  In  this  letter,  the  density  functional 
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theory  (DFT)  calculations  are  carried  out  by  DMoP  package  [22],  The  DFT  calculations  are  performed  by 
using  generalized  gradient  approximation  (GGA)  with  the  functional  parameterized  by  Perdew,  Burke,  and 
Enzerhof  (PBE)  [22],  And  the  double-numerical-polarization  (DNP)  basis  set  is  employed.  Self-consistent 
field  calculations  are  done  with  a  convergence  criterion  of  10“®  hartree  on  the  total  energy.  The  Brillouin 
zone  integration  is  performed  using  a  6x6x6  Monkhorst-Pack  k-point  grid.  All  the  structures  are  fully 
optimized  with  a  convergence  criterion  of  0.002  hartree/A  for  the  forces  and  0.005  A  for  the  displacement. 
A  real-space  cutoff  of  4.0  A  for  the  atom-centered  basis  set  is  chosen  to  increase  computational  efficiency 
while  not  significantly  affecting  the  magnitude  of  inter-atomic  forces  or  the  total  energies.  The  Gaussian 
smearing  of  electron  density  is  applied  with  the  energy  range  of  0.1  eV. 

The  electrical  conductivity  a,  the  thermal  conductivity  due  to  electrons  Xg,  and  the  Seebeck  coefficient  S, 
are  obtained  from  the  electronic  structure  with  the  solution  of  1 -Dimensional  Boltzmann  transport  equation 
as  mentioned  above.  Here  we  use  the  experimental  data  of  n-type  alloy 

{n  =  1.5x10^“ cm , cr  =  3400(5'cm“')  from  Ref.  [14],  and  the  dependence  relation  [19]  between  the  mobility 
and  carrier  concentration  in  alloys  to  calculate  the  carrier  concentration  dependent  relaxation  time. 
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5.  Results  and  discussions 


5.1  Thermal  electric  properties  of  isotope  doped  silicon  nanowires 

In  this  part  we  present  our  results  on  size  dependence  of  physical  properties  of  Si  NW. 

Fig.  1(a)  and  (b)  show  the  size  effect  on  o  and  S  with  different  electron  concentration,  o  increases 
slightly  as  the  diameter  increases,  while  the  Seebeck  coefficient  S  decreases  remarkably.  The  size 
dependence  arises  from  quantum  confinement  effect  on  the  electronic  band  structure. 

Figure  1(c)  and  (d)  show  that  S  decreases  as  the  carrier  concentration  increases,  while  o  increases 
as  more  carriers  are  available  to  transport  charge.  Detailed  discussions  can  be  found  in  our  publication  in 
APL  (2009). 

The  power  factor  P  (P  =  S^cr),  is  an  important  factor  influencing  the  thermoelectric  performance 
directly.  Figure  2  (a)  shows  the  power  factor  versus  carrier  concentration  there  is  an  optimal  carrier 
concentration  A/^ax  yielding  the  maximum  attainable  value  of  Pmsx 

The  figure  of  merit,  ZTis  another  important  thermoelectric  characteristic.  In  the  calculation  of  ZT, 
both  electron  and  phonon  contribute  to  the  total  thermal  conductivity.  Figure  3(a)  shows  the  dependence  of 
electron  thermal  conductivity  on  carrier  concentration  for  different  transverse  size.  It  is  clear  that  electron 
thermal  conductivity  increases  with  carrier  concentration  as  it  is  proportional  to  electronic  conductivity,  and 
SiNW  with  larger  diameter  is  with  higher  electron  thermal  conductivity.  For  SiNWs  with  length  in  pm  scale, 
the  phonon  thermal  conductivity  increases  with  diameter  increases  remarkably  until  the  diameter  is  larger 
than  about  hundreds  nms. 
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Figure  1.  (a)  Electrical  conductivity  vs  cross  sectional  area  with  different  carrier  concentration,  (b)  S 
vs  cross  sectional  area  with  different  carrier  concentration,  (c)  Electrical  conductivity  vs  carrier 
concentration  with  fixed  cross  section  area  of  1 .1  nm^.  (d)  S  vs  carrier  concentration  with  fixed  cross 
section  area  of  1 .1  nm^. 


Figure  2.  (a)  Thermal  power  factor  of  SiNW  vs  carrier  concentration  with  three  different  transverse 
dimensions,  (b)  Maximum  power  factor  vs  cross  sectional  area,  (c)  Nmbx  vs  cross  sectional  area,  (d) 
Size  dependence  of  the  maximum  room  temperature  cooling  power  density  of  SiNW  with  length  of 
1pm. 
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Figure  3.  (a)  Thermal  conductivity  due  to  electrons  vs  carrier  concentration  for  SiNWs  with  different 
transverse  dimensions,  (b)  ZT  vs  carrier  concentration  for  different  isotope-doped  SiNWs  (^®Sii. 
x^^Six  NWs)  with  fixed  cross  section  area  of  2.3  nm^.  (c)  ZTuax  vs  the  concentration  of  ^®Si  atom,  (d) 
AZ/wax  vs  the  concentration  of  ^®Si  doping  atom 


5.2  Thermal  electric  properties  of  Ge  doped  silicon  nanowires 

In  Figure  4(a)-(d)  we  show  the  Ge  content  dependence  of  a  and  S  with  different  electron 
concentration  for  NWs.  Here  we  choose  the  heavily  doped  case  because  the  maximum  ZT  is 

expected  to  occur  for  high  carrier  concentration,  which  has  been  confirmed  experimentally,  and  also  to 
reduce  the  unknown  effects  of  the  band-gap  uncertainty  on  the  calculated  values  of  S. 

Figure  5  (a)  and  Figure  5(b)  show  the  power  factor  P  {P  =  xa)  versus  carrier  concentration  for 
n-type  and  p-type  wires.  There  is  an  optimal  carrier  concentration  Nmbx  yielding  the  maximum  attainable 
value  of  P/wax-  Here,  W/wax  is  around  «  =  2.0 x  10^“ for  n-type  wires  and  «  =  1.0 xl0^'cm"^for  p-type 
wires.  As  the  Ge  content  increases,  for  n-type  wires,  P/wax  decreases  remarkably  until  x=0.5,  and  then 
increases  with  x  as  shown  in  Figure  5(c).  For  p-type  wires,  P^ax  decreases  (Figure  5(d))  with  the  Ge 
content  increases  monotonically  until  x  =  0.8 .  The  large  fluctuation  is  partially  due  to  the  effects  of  the 
band-gap  uncertainty  on  the  calculated  values. 

The  most  significant  results  in  the  current  work  are  shown  in  Figure  6.  Here  we  focus  on  the  relative 
change  of  ZTin  n-type  NWs  and  compare  it  with  pure  Si  NWs.  Using  our  calculated  power 

factor  from  Figure.  5,  and  the  relative  phonon  thermal  conductivity  ( iCsiOe^^si  -  here  and  /f^^.g^are 
thermal  conductivities  of  SiNW  and  Si^_^Ge^  NW,  respectively)  calculated  by  using  molecular  dynamics 
method  from  Ref.  1 1 ,  the  dependence  of  ZT^.^  on  Ge  content  x  is  shown  in  figure  4.  The 

ZTsi^  Ge  values  increase  with  Ge  content,  reach  a  maximum  and  then  decreases.  At  low  Ge 
content,  the  small  ratio  of  Ge  atoms  can  induce  large  increase  in  ZT.  For  instance,  in  the  case  of 
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Sio.8Geo.2  NW,  namely,  20%  Ge,  ZT^.^  / ZT^^  is  about  3.  And  with  50%  Ge  atoms  (Sio.sGeo.s  NW),  the 

ZTsi^  Ge  !  Si  high  as  4.3.  The  similar  dependence  of  ZT  is  also  observed  in  p-type  wire, 

although  its  ZT  is  much  lower  than  that  in  n-type  wire.  This  phenomenon  can  be  understood  as  the 
following.  ZTis  contributed  by  both  power  factor  and  thermal  conductivity. 


Figure  4.  a  vs  Ge  content  x  with  different  carrier  concentration  for  n-type  wires  (a)  and  p-type  wires  (b).  S  vs 
Ge  content  with  different  carrier  concentration  for  n-type  wires  (c)  and  p-type  wires  (d). 
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Figure  5.  Thermal  power  factors  of  Si^_^Ge^  NWs  versus  carrier  concentration  with  three  different  Ge 

contents  for  n-type  wires  (a)  and  p-type  wires  (b).  Maximum  power  factors  versus  Ge  content  for  n-type 
wires  (c)  and  p-type  wires  (d). 


Figure  6.  Z4 
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